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ABSTRACT 

We perform radiation hydrodynamics simulations on the evolution of galactic gas disks irradiated by 
ultraviolet radiation background. We find gas disks with Ah ^ lO^-'^cm"^ exposed to ultraviolet radiation 
at a level of I21 — 1 can be self-shielded from photoheating, whereas the disk with Ah < lO^-'^cm"^ 
cannot. We also find that the unshielded disks keep smooth density distribution without any sign of 
fragmentation, while the self-shielded disks easily fragment into small pieces by self-gravity, possibly 
followed by star formation. The suppression of star formation in unshielded disks is different from 
photoevaporation effect, since the assumed dark halo potential is deep enough to keep the photoheated 
gas. Presence of such critical threshold column density would be one of the reason for the so-called 
down-sizing feature of present-day galaxies. 

Subject headings: galaxies: formation — radiative transfer — hydrodynamics 



1. INTRODUCTION 



Recent advances of available computational resources 
enable us to simulate the entire disk galaxy to investi 
gate parsec scale structures in the diskdWada fc Norma: 
2007t iRobertson fc Kravtsovll2007t iTasker fc Brvanll200 



Saitoh et al.l 120081 ). In addition, numerical techniques on 



Radiation HydroDynamics (RHD) simulations are now 
being developed, e specially in the field related to the 
cosmic rcionizationtone diii fc Abel Il200]| _ Mellema et al.l 
2006; Riikhorst et al. 2 00(1 ISusa Il200(il: llliev et all 120061: 
^shida ct al. 2007; jQiu et all 120071 : IWhalen fc NormanI 
2008; Altav ct al. 20081). Using these schemes, we are now 



ready to tackle the numerical simulations of the detailed 
structures in galactic disks coupled with radiation transfer. 

The importance of ultraviolet radiation transfer ef- 
fects on the galaxy formation/star formation in galaxies 
have been pointed out by several authors. First of all, 
the observed star formation threshold in galactic disks 
(|Kennicuttlll989l; iMartin fc Kenn icu"ttll200lh could be ex- 
plained by the self-shielding effects of galactic dis k s from 
externa l ultr aviolet radiation field ()Schavel 1200 iL l2004f) . 
ISchavd ()2004[ ) obtained the physical state of the disks ex- 
posed to external ultraviol et field, using one dimensional 
calculation with CLOUDYfFerland"2000'). Based upon the 
liner stability arguments (Toomre 1964 ), he has demon- 
strated that the inner parts of the galactic disks shielded 
against external radiation field by dust absorption could 
be gravitationally unstable, whereas the unshielded outer 
parts of the disks are stable. These results nicely explain 
the observed features of galactic disks, however, it has not 
been tested by multi-dimensional RHD simulations. 

Secondly, in the context of galaxy formation, radia- 
tive feedback by ultraviolet background is expected to 
be very important especially for low mass galaxies. The 
photoheating can heat the gas up to lO^K, which pre- 
vent the gas from collapsing, in case the gravitational 
potential of the dark matter halo is not deep enough 



to retain the photoheated gas. Such feedback effects 
are quoted as ph otoevaporation, which is well stud- 
ied at various levels fUmemura fc Ikeuchi"l984';'Efstathiou| 
1992; Babul fc Rces 1992; Thoul fc Weinberg 1996- 
Barkana fc Loe b 1999; iFerrara fc Tolsto v 2000; Gnedinl 
20001 iKitavama et al.l l2000l 120011: ISusa fc Umemm^ 
2004a|P~~ In addition, there are some evidences in 
dwarf galaxies that the star formation is suppressed in 
"hot" (~ lO^K) phase (Young fc Lo 1997a, b), which in- 
fer that photoheating can suppress the star formation in 
galaxies, even if the potential of the dark halo is deep 
enough to prevent the gas from evaporating. This issue 
could have a concern to the so-calle d down-sizi ng prob- 
lem i n nearby galaxies (|Cowie et al.ll9 96: Kauffm ann et al.l 
|2003| ). They found that old galaxies are massive, while 
the young ones are less ma ssive. Such trend continues to 
higher redshiftfe.g ., Koda ma et alll2004l ; lE^b et al.l [20061 ; 



iReddv et al1 l2006; Papovic h et al.H2006f) . One of the pos- 
sible interpretation of this feature is that star formation 
proceed very rapidly in massive galaxies, whereas it is 
a slow process in less massive galaxies for some reason. 
I Kauffmann et al] (j2003) also pointed out that these two 
groups are well defined, i.e. they have a clear boundary at 
^ 10^°Mq in stellar mass. This critical mass scale is too 
large to be related to the photoevaporation mechanism, 
however, photoheating could still be the candidate to ex- 
plain the down-sizing mass if we can show the mechanism 
to suppress the star formation even in halos in which the 
photoionized gas can be inherent. 

In this paper, we examine the fragmentation of gas disks 
embedded in dark halo potential under the ultraviolet 
radiati on field, us ing the recently developed RHD code 
RSPH (|Susa II2006I ). This paper is organized as follows. In 
the next section, the numerical scheme is briefly summa- 
rized. In section [21 the setup of numerical simulations are 
described. We show the results of numerical simulations 
as well as the analytic estimate in section 21 Sections [5] 
and [6] are devoted to discussions and conclusion. 
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2. METHODOLOGY 

We perform n umerical sim ulations by Radiation- 
SPH(RSPH) code (|Susa Il2006f) . The code can compute 
the fraction of primordial chemical species e~, H^, H, H~, 
H2, and by fully implicit time integration. It also can 
deal with multiple sources of ionizing radiation, as well as 
the radiation at Lyman- Werner band. 

Hydrodynamics is calculated by Smoothed Particle Hy- 
drodynamics (SPH) method. We use the version of SPH by 
[Um cmura (1993) with the modification on the SPH ker- 
nel function and symmetrization of equ ation of motions 
according to ISteinmetz fc Miilleil ( 19931). We also adopt 
the particle resizing formalism by iThacker et al. I (|2000[ ) 
in which the number of neighbor SPH particles are kept 
almost constant without sudden changes. 

The non-equilibrium chemistry and radiative cooling for 
primordial gas are calcu lated by the code developed by 
ISusa fc Kitavamal (j2000[ ). where H2 cooling and reaction 
rates are mostly taken from lGalli fc Fallal ( 19981 ). As for 
the photoionization process, we e mploy so-called on-the- 
spot approximation ()Spitzeijri978l ). In this paper, we also 
added the radiati ve cooling due to metals em ploying the 
formula given in Dalgarno fc McCravl (|l972f). assummg 



Z — 10 '^Zq. The metallicity of nearby disk galaxies are 
normally larger than Z = 10~^Zq, however, the observed 
metaUicity of t he IGM is Z - IO^^Zq - IQ-^Zq (e.g., 
iSongailal |2001[ ) which is the ingredient of the galaxies. 
Since we are interested in the forming disk galaxies, we 
employ Z = 10~'^Zq in the present set of simulations. Re- 
mark that the radiative transitions of heavy elements are 
not the dominant process of r adiative cooling for > lOO OK 
(|Boehringer fc Hensleil [19891 ISusa fc Umemurall2004al) in 
case Z = 10~'^Zq is assumed. H and H2 molecules are the 
main coolant of such gafl. 

The optical depth is integrated utilizing the ncigh- 
bor lists of SPH par ticles. In our old scheme 
(|Susa fc Umemu ra''2004a'). we create many grid points on 
the light ray between the radiation source to an SPH par- 
ticle. In the present scheme, we do not create so many 
grids. Instead, we create one grid point on the light ray 
per one SPH particle on its upwind. We find the neighbor- 
ing 'upstream' particle for each SPH particle on its line of 
sight to the source, which corresponds to the grid point. 
Then the optical depth from the source to the SPH par- 
ticle is obtained by summing up the optical depth at the 
'upstream' particle and the differential optical depth be- 
tw een the two particles. The more details are described 
m ISusa I (I2006D . We assess the optical depth for ionizing 
photons as well as the Lyman- Werner band photons by 
the method described above. In the present version of the 
code, the dust opacity is not included in the calculation. 

The c ode is tes ted for various standard RHD 
problems ljSusa I [2006t ). and already applied to the issues 
on the radiative feedback effects of first generation stars 
(|Susa fc Umemurall2006l : ISusall2007[ ). We also take part in 
the code com parison pro ject w ith other radiation hydrody- 
namics codes ()Iliev et al.ll2006f ) in which we find reasonable 
agreements with each other. 

3. SETUP OF NUMERICAL SIMULATIONS 



We perform numerical simulations of galactic gas disk 
embedded in a dark halo potential. The dark halo poten- 
tial is fixed as 
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= V2 — lOOkm/s. 



where r denotes the distance form 
Ikpc, ai = 0.3kpc, 02 = 5kpc and vi 
Th us the assumed po t ential is similar to the one employed 
by IWada fc NormanI (|2007l ). except that the rotation ve- 
locity at given radius is smaller by a factor of 2, since we 
are interested in the forming galaxies relatively less mas- 
sive than our Galaxy. 

In the present set of simulations, we assume initially 
uniform disk (Fig. [T]) with slight perturbations (displace- 
ments of SPH particles are < 10 %). The initial thickness 
and radius of the disks. Hi, i?disk are lOOpc and 3kpc for 
all of the runs, respectively. The initial velocity of the gas 
particles in the disk are assumed to be the Kepler rota- 
tion velocity around the dark halo potential center. Using 
equation ([T|) , the rotation velocity at the edge of the disk 
r — i?disk is as large as ~ lOOkm/s. Thus, the dark halo 
potential is deep enough to retain the gas even if it is pho- 
toheated to ^ lO^'K. The initial densities of the disks, pi, 
are changed depending on the models, those are listed in 
Tabled 

The light rays of external ultraviolet radiation field are 
assumed to be perpendicular to the gas disk. We consider 
two directions (upward/downward) as shown in Fig. [TJ 
This particular choice of the ray-direction is an approxi- 
mation if we consider background radiation field exactly, 
since background photons are coming from all angles. If 
we take into account such effects, the boundary between 
the shielded and unshielded regions are softened, which 
might change the self-shielding effects slightly. However, 
as far as we keep the number of photons coming into the 
disk unchanged, the self-shielding condition is also almost 
unchanged, since it is determined by the balance between 
the number of recombining hydrogen atoms and the num- 
ber of photons. The flux of the field at the Lyman limit 
is given so that the mean intensity at the midplane of the 
disk equals to I21 x 10~^^ergs~^cm^2Hz^^str^^ in case the 
opacity of gas disk is ignored. We consider an AGN type 
spectrum, which is proportional to v^^. 

For all of the runs, we employ common physi- 
cal/numerical parameters of gas disk listed in Tabl^ 
while the model dependent parameters are listed in Table 
[2I As shown in the table, we perform ten runs changing 
the density of the disk and the intensity of the radiation 
field. 

Since the resolution of the simulation is constrained by 
the mass of an SPH particle, we have the critical density 
above which the gravitational fragmenta tion are not cap- 
tured properly. According to B ate fc BurkertI (|1997D . the 
density nH,ros is given as 
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where mgph, Mici, ?inin denote the mass of each SPH par- 
ticle, number of neighbor particles and the minimal tem- 



^ We have to keep in mind that for low temperature (T < 500K), or high density (njj ^ 10*cm 
to decide the fate of cooling gas (.Omukaii, 2000,) . 



realm, even Z ~ 10 ^.^0 could be important 
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perature set in the simulations, respectively. The sym- 
bols /cb,G, rrip have ordinary meanings, i.e., they repre- 
sent Boltzmann constant, gravitational constant and pro- 
ton mass, respectively. In the present set of simulations, 
"^H.ios equals to 235 cm~^, except the runs B/2,Br/2,B/8 
and Br/8 in which the mass of a SPH particles is larger 
than other regular runs by a factor of two or eight (see 
Table [2]). The softening length of gravitational interaction 
between SPH particles , e, is set so as to satisfy 

^ _ 1 / SiV^loimSPH \ .gS 

2 v47rmp7iH,rcs/ 
This expression guarantees that the self-gravity of a dense 
clump above rtH.rcs is softened. 

Given these configurations, we can assess the Toomre's 
Q value to understand the stability of the initially uniform 
disk. Toomre's Q is defined as 



Here Cs,k and S denote the sound speed of the gas, 
epicyclic frequency and the surface density of the disk. 
In case Q > 1 holds, the disk is gravitationally stable be- 
cause of the Coriolis force/t hermal pressur e, whereas it is 
unstable if Q < 1 is satisfied (jToomrelll964[ ). 

The actual Q parameters for three models A, B and C 
(see Table [2|) are plotted in FigO Here Cs is assumed to 
be the sound speed corresponding to Tmin = 300K. Since 
we assume uniform disks for all models, the disks are more 
unstable at the outer radii. In the present simulations, the 
radius of the disks are 3kpc for all models. Thus, there ex- 
ist critical radii above which the disks are unstable if those 
are cooled down to Tj^in, although the unstable region is 
narrow for model A ( it is unstable only at the edge of the 
uniform disk). It is worth noting that if we employ Tmin 
lower than 300K, the unstable region should expand, since 
the gas temperature will go down to Tmin, that results in 
smaller Q. 

4. RESULTS 

4.1. Stability of the disks without Radiation 

First, we present the results for /21 = 0, i.e. no ex- 
ternal radiation field. The face-on/edge-on view of the 
snapshots for the disks in models A, B and C are shown in 
FigO The snapshots are taken at i = SOOMyr for model 
A, t = 120Myr for model B and t = 40Myr for model C. 
The particular choices of these times when the snapshots 
are taken basically corresponds to the time when the time 
steps at dense gas clumps collapsed below nn.rcs becomes 
so short ( < lO^yr ) that physical time in numerical com- 
putation evolve very slowly. Since the growth time scale of 
gravitational instability is shorter in denser disks, we take 
the earlier snapshot in model C than others. 

It is clear that the disks in all models fragment into 
small filaments/knots. At the same time, the inner parts 
of the disks are stable. As shown in Fig [21 Q values are 
larger than unity at smaller radii, whereas they are smaller 
than unity at larger radii. Therefore, the inner parts of the 
disks are stable, which is consistent with the present re- 
sults. In addition, the critical radii outside which the disks 
are unstable predicted by Toomre criterion are also con- 
sistent with numerical results. Indeed, the critical radii 



read from FigOare 2.8kpc, 1.6kpc and 0.8kpc for models 
A,B and C, those are almost consistent with the boundary 
radii of smooth parts of the disks. 

In Figdl the phase diagram (nn — T plane) of the three 
runs at the same snapshots as in Fig|3] are shown. The 
gas temperature in all of the models are cooled efficiently, 
which are almost close to Tmin = 300K at high density 
realm (nn ^ lO'^cm"'^). Such efhcient cooling justify the 
stability argument based upon the Q values plotted in 
Fig (21 in which T = T,„in is assumed. 

These phase digram also show that very dense regions 
above nn.rcs sue formed in all of the models. These dense 
regions correspond to the fragments presented in Fig. [31 

4.2. Effects of Ultraviolet Radiation on the Disk Stability 

Now, we show the snapshots of the runs with I21 = 1 
(Ar, Br and Cr) in FigO i.e. t = SOOMyr for model Ar, 
t = 120Myr for model Br and t = 40Myr for model Cr. 
The snapshots are synchronized to the corresponding ones 
in FigO At a first glance we find clear difference from 
the runs without ultraviolet radiation. In model Cr, we 
find small filaments and knots as were found in model C, 
whereas we cannot find any sign of fragmentation in mod- 
els Ar and Br. In addition, the density of the disks in 
models Ar and Br becomes lower than that in models A 
and B. We also observe that the disks are geometrically 
thicker than the previous ones. In Fig|6l the phase dia- 
grams are shown for models Ar, Br and Cr. In all of the 
runs, significant amount of materials are photoheated up 
to > lO^K in low density region, that contribute to the 
disk thickening. Additionally, no very dense gas compo- 
nent above riH,rcs is found in models Ar and Br, while we 
find it again in model Cr as was found in model C. Thus, 
these results indicate that the presence of external ultravi- 
olet radiation field suppress the fragmentation of the disks 
with low density models (Ar and Br), although it does not 
for high density model (Cr). 

4.3. Self- Shielding from Ultraviolet Radiation 

The difference between low density models(Ar, Br) and 
high density model(Cr) comes from the self-shielding ef- 
fects. FigjTl shows the density (top), temperature (mid- 
dle) and HI fraction j/hi (bottom) distributions along the 
direction perpendicular to the disks. In model Cr (right 
column), the gas near the midplane (z = 0) is neutral, 
cold and dense, because of the self-shielding. On the con- 
trary, in models Ar and Br, cold and dense regions near the 
midplane is relatively smaller than that in model Cr. To 
be more quantitative, FigjSl shows the temperature prob- 
ability distribution functions (PDFs)for models Ar,Br,Cr 
(top) and A,B,C(bottom). The temperature PDFs show 
the mass fraction of the disk found within a logarithmic 
temperature bin A(logT) — 0.06. In the runs without ra- 
diation field (A,B and C), most of the mass is condensed 
in the coldest phase at Tmin- On the other hand, the tem- 
perature PDFs are broader in the runs with radiation. In 
models Ar and Br, most of the mass in the disks exist in 
the hot phase {'^ 10'' K), while cold phase (< lO'^K) is 
dominant in model Cr. Such large difference in tempera- 
ture distribution results in different stability of the disks. 
According to Fig. |21 Q values for models Ar and Br satisfy 
> 0.5 if r < 3kpc and T = Tmin- In actual simulations, we 
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find T ^ lO^K, which means that Q values are larger by a 
factor of yiWsOO = 5.77 than 0.5, i.e. Q > 1 is achieved 
everywhere in the disks for these models. Thus, these disks 
are stable, while the disk in model Cr is unstable because 
of its coldness. 

A rough estimate of self-shielding criterion for primor- 
dial gas has been derived by Susa fc Umemural (|2000af) . 
who evaluate the shielded photoheating rate to compare 
with the peak H2 cooling rate below lO'^K. Remark that 
heavy elements are not the dor ninant coolant for lO^K < 
T < lO^K, in case Z < 10-2^0 dSiisa & Umemurall2004a^ 

Similar argument is also found in Corbelli e t al. (1997), 
which describes th e thermal instabili t y of th e primordial 
gas. According to ISusa &: Umemural (|2000aD . the photo- 
heating rate per unit volume at the midplane of the disk 
for the optically thick limit is given as ?T.H.cyHi,cW where 



3 l + /3^"^ 



(5) 



Here nn.c and ym,c are the number density of the hydrogen 
nucleus and HI fraction at the midplane of the disk, a^,^ is 
the photoionization cross section at Lyman limit, z^l is the 
Lyman limit frequency, denotes the incident ultraviolet 
intensity at Lyman limit, T{x) is the gamma function. /3 
used in this equation is defined as /? = l+{a — 1) /3, where 
a denotes the spectral index, i.e., I^, oc j/^" is assumed. 
Tj^L denotes the optical depth at Lyman limit, which is 



written as. 



(yni) Nna^ 



(6) 



Here denotes the column density of the disk, (?/hi) 
is the HI fraction averaged along the direction per- 
pendicular to the disk. The peak H2 cooling rate 
per unit volume below 10*K is described as ^Ch2 , 
where Ch-, ^ lO'^'^erg s~^ cm^ ()Shapiro fc KangI il987t 
ISusa fc Umemur£ill2004a[) . In addition, we assume hydro- 
static equilibrium in the vertical direction. As a result, 
the number density of hydrogen nucleus at the midplane 
of the disk is related to the column density as follows: 



nGN^fimp 
2^! ■ 



(7) 



Here the gas disk is assumed to be isothermal with sound 
velocity Cg. Letting the photoheating rate be equal to the 
H2 peak cooling rate, combined with the equations ([S]), © 
and ([7|), we obtain the threshold column density iVH,sh, 
above which the photoheating is shielded enough for the 
midplane to cool down to < lOOOK: 

\iGpLmpCu2ym.cKVm) (/3+1)/ 
Using the present assumptions (a = 1, /21 = 1), we have 



TVH.sh =i 2.5 X lO^^cm"^ 
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Here we also assume j/hi.c = 1 hi the above expression, 
since it is true for all models Ar, Br and Cr. 



On the other hand, the initial column density of the 
disks in models Ar, Br and Cr are 6.3 x 10^°cm~^, 
1.3x lO^^cm^-^ and 3.8 X lO^^cm"^, respectively. Thus, the 
numerical simulations indicate 1.3 x lO^^cm^^ < iVn.sh ^ 
3.8 X lO'^^cm"^, which almost agrees with the rough ana- 
lytic estimate shown above. 

4.4. Probability Distribution Functions of Density Field 

The differences among the six models (A,Ar,B,Br,C and 
Cr) are clarified in Fig. [9] in terms of density proba- 
bility distribution functions (PDFs). Six panels represent 
the snapshots of the density PDFs for six models. Here 
the density PDFs are defined as the mass fraction found 
within a logarithmic density bin A(lognH) = 0.18. The 
histograms drawn by thin lines represent the initial PDFs, 
while the thick lines show the evolved ones. 

It is clear that in the models A,B,C, and Cr, the PDFs 
extend beyond the resolution limit nn.res- Thus, dense self- 
gravitating fragments formed in these disks, which would 
lead to the star formation activities. 

It looks strange at first that very dense clumps form 
beyond nn.rcs, although the gravitational force is softened 
for riH > J^H.res (sec eq. ([3])). In fact, even if gravitational 
force could be neglected at nn ^ ?T-H,res, matters accrete 
from the outer envelope. As a result, density (or pres- 
sure) of the clump becomes higher than nH,rcs, in order 
to compensate for the increasing ram/thermal pressure of 
accreting material. 

On the other hand, in the models Ar and Br, PDFs 
have very sharp cut-off below nH,rcs for the snapshot at 
t = 300Myr(Ar) and t = 120Myr(Br), although models 
A and B have clear fragmentation signature at the same 
snapshots. Thus, the impression from the montage of the 
disk is correct, i.e. the disk in models Ar and Br do not 
fragment into dense clouds. This result could be inter- 
preted that the star formation in these disks are heavily 
suppressed by ultraviolet radiation field. 

4.5. Convergence of PDF 

The convergence of density PDF is checked for models 
B and Br. We perform runs with larger SPH particles 
mass by a factor of two and a factor of eight. These runs 
are tagged as B/2, Br/2 and B/8, Br/8 (Tabled. FigdH] 
shows the density PDFs for models B, B/2, B/8 (top) and 
models Br,Br/2,Br/8 (bottom) at t = 40, 100, 120Myr. In 
the half resolution runs (B/2 and Br/2), the density res- 
olution limits are four times smaller than those in B and 
Br, because nH,ros inversely proportional to the square of 
mass resolution (see equation Similarly, the density 

resolution limits in B/8 and Br/8 are 64 times smaller than 
B and Br. The PDF histograms basically agree very well 
below nH,rcs as expected. 

In the runs B/8 and Br/8, the density resolution limits 
are so low that the high density regions (> 100cm~^) are 
not captured properly. As a result, it seems to be difficult 
to distinguish the PDFs for models B/8 and Br/8 from 
each other. On the other hand, the results from the runs 
B/2 and Br/2 indicate the same conclusion as we found in 
the runs B and Br, since the density resolution is enough 
to capture the self-shielded dense regions. Thus, the disk 



The expression "H is 2 times larger than the equation (A12) in ISusa fc Umemural (I2000al) . since the radiation field irradiate the right side of 
the disk as well as the reverse side in the present simulations. 
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in B/2 fragments due to gravitational instability, whereas 
that in Br/2 is stable. Therefore we conclude that the 
resolution of the present regular simulations are enough 
to capture the fragmentation of the disks under the as- 
sumptions we employed, whereas in lower resolution runs 
such as B/8 and Br/8 are not able to describe the present 
physical situation. 

5. DISCUSSION 

It should be emphasized that present results di- 
rectly prove the suppression of star formation ac- 
tivities by ultraviolet background radiation in halos 
with Tvir > 10'* if. It has already been known that 
the star formation is suppressed in less massive ha- 
los with Tvir ^ 10'*K, since the gas in the dark 
halo potential evaporates d ue to the high thermal 
pressure o f photoheated gas fUmemura & Ikeuchi' '1984, 
Efstathioul ^1992: Babul & Rccs 1992; Thoul & Weinberg 
1996; Barkana & Locb 1999|; IFerrara fc Tolstovl |2000[ 
Gnedin ,,2000.) if the s e lf-shielding effect is not i mport ant 
( Kitavama etalll200d 120011 : ISusa k Umemurall2004al lbl). 
On the other hand, in the present simulations, the gas in 
the dark halo do not evaporate because of the deep grav- 
itational potential. We find that even in such halos, star 
formation could be heavily suppressed in case gas is con- 
figured to form disks with low column density. 

Another important result found in the present calcula- 
tion is the presence of a clear boundary in the column 
density of the disk, below which star formation is heav- 
ily suppressed. According to the numerical results, the 
critical column density is TVh ^ 1 — 4 x lO^^cm"^. This 
threshold is very interesting mainly for two reasons. First, 
although the present simulations are performed with fixed 
dark halo mass, we can try to convert the critical column 
density into dark halo mass. The mass of the uniform disk 
with given TVh and disk radius i?disk is 

Therefore, if the mass of the host dark halo is 7(~ ^I-q/^) 
times the disk mass, the critical dark halo mass is ^ 
4 X 1O^"A/0. This value is still s maller than the critical 
scale found by iKauffmann et al.l (|2003D by a few factor, 
but if we take into account the internal feedback effects 
such as UV radiation from internal sources (AGN, mas- 
sive stars) or supernova, the threshold might account for 
the observed critical down-sizing mass. We also point out 
that if we start the simulation from cosmological setup, the 
threshold mass could be raised more, since photoheating 
might be able to penetrate deeper into the disk because of 
initially less dense configurations. Thus, inclusion of such 
physics as well as starting simulations from cosmological 
setup are necessary to obtain more precise understanding 
of star formation in forming disk galaxies. 

Secondly, the critical column density is as large as 
the star formation t hreshold found in local d isk galax- 
ies (|Kennicuttl [19891: iMartin fc Kennicutil200lh . In fact, 
the observed star for mation t hreshold of galactic disks has 
been investigated bv lSchavi ()2001l . [200^ from theoretical 
side. He suggest that the galactic disk could be gravita- 
tionally stable if it is unshielded from the external UV 
radiation, since Q value of the photoheated disk easily 



exceeds unity. As a result, gravitational fragmentation 
of the di sk is suppressed, so the star formation activi- 
ties do. ISchave fc Dalla Vecchial (|2008D performed sim- 
ple hydrodynamics simulations taking into account the 
star formation threshold column density using effective 
equation of state in multiphase medium, although they 
do not solve radiation transfer equations explicitly. They 
found the threshold column dens ity is as l arge a s 4M0pc~^ 
(Fig.4 and 5 in Schave fc Dalla Vecchial (|2008D ). which is 
5 X 10^°cm~^ in cgs unit. This value is smaller than that 
obtained in this paper by a few factor('^ 2.5 x lO^-'^cm^). 
The basic mechanism they propose to suppress the star 
formation of the disk is same as the one found in the 
present simulations, except that we do not include the ef- 
fects of dust extinction required especially for present-day 
disk galaxies. The gas disk is more easily self-shielded by 
dust absorption, since the dust opacity at the Lyman limit 
frequency is as large as HI continuum for solar metallic- 
ity. Thus, it is it is reasonable th at we have larger critical 
column density than obtained in ISchave fc Dalla Vecchial 
(|2008f ). On the other hand, in the present simulations, 
we assume Z = 10~^Zq. Therefore, the dust opacity at 
Lyman limit is smaller than HI continumm opacity by 
two orders of magnitude. As a result, dust opacity has 
much smaller impact on the self-shielding effects at such 
low metallicity. In any case, we would dare to mention 
that present results almost succeeded to probe the pres- 
en ce of sta r form ation threshold column density proposed 
bv ISchavi (|2004l ) utilizing the full 3D radiation hydrody- 
namics simulations. More realistic calculations for present- 
day galaxies including the effects of dust extinction are left 
for future works. 

It is also worth noting that star formation rate (SFR) of 
the disk above the threshold column density is consistent 
with the observed value, although it is difficult to compare 
the present results directly with observed SFR, since we 
do not take into account the local stellar feedbacks, the 
effects of dust particles, metal enrichment followed by the 
radiative cooling by abundant metals. Despite of such is- 
sues, SFR cloud be obtained by the density PDF assuming 
the stars are formed in dense clumps within local free-fall 
time. It is given as 

S* = Cca/ G/3cSdisk/c (11) 

where is the star formation efficiency, Pc is the threshold 
density above which the gas is converted to stars, Sdisk 
denotes the surface density of the disk, and fc denotes 
the mass fraction of gas in the disk condensed into dense 
clumps with p > pc- /c is the quantity which can be ob- 
tained by integrating the density PDF above pc- In the 
present set of simulations with radiative feedback, model 
Cr is the only one in which star formation is expected 
since self-shielded cold fragments emerges in the run. If 
we use e = 0.1 and nH.ros to assess the threshold den- 
sity Pc, we obtain fc ~ 0.066. Using the surface den- 
sity of the disk in model Cr, the SFR is evaluated as 
S* — O.O3M0yr~-^kpc~^. Although we cannot discuss the 
dependence of S FR on Sdisk, the value is consistent with 
the observation f Kennicuttlll998l ). This reasonable agree- 
ment infer the validity of the present numerical models. 

We also point out an o ther s tandpoint on this is- 
sue by ISusa fc Umemural (|2000af ). in which they sug- 
gested that self-shielding from ultraviolet background 
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could be a key mechanism to deter mine the morphology 
of galaxies ljSusa fc Umemurall2000bD . although their argu- 
ments based upon 1-dimensional radiation hydrodynamics 
calculations. Unfortunately, it is impossible to relate the 
present results to morphology bifurcation of galaxies, since 
we assume a disk by hand in our simulations performed so 
far. From this point of view, again we need to perform 
simulations from cosmological initial conditions. 

6. CONCLUSION 

In this paper, we perform radiation hydrodynamics sim- 
ulations on the fragmentation of galactic disks under the 
ultraviolet radiation background. We find that ultraviolet 
radiation field strongly suppress the star formation in the 
disks in case the photoheating is not shielded enough. We 
emphasize that this suppression is different from photo- 
evaporation effect, because the rotation velocities at the 
outer boundary of the disks in the present set of simula- 
tions are ^ lOOkm/s, which are fast enough to retain the 



photoheated gas. In our simulations, we find a threshold 
column density of the disk (~ lO^^cm"^) above which the 
fragmentation is not suppressed. It is similar to the star 
formation threshold column density observed in nearby 
galaxies. Presence of such critical threshold would be 
one of the reason for the so-called down-sizing problem 
in nearby galaxies. 
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Table 1 
Common Parameters 



symbol 


numerical value employed 


description 


N ■ 

^ ^nei 


50 


Number of neighbor particles 


Hi 


lOOpc 


Initial disk thickness 


-Rdisk 


3kpc 


Initial disk radius 




lO^K 


Initial disk temperature 


T ■ 

min 


300K 


minimum temperature 



Table 2 

Model dependent Parameters 



Label 


hi 


Pi[M0/pc-^ 


simulated time 


# of SPH particles 


'^H.res 


cm 


e[pc 


A 





0.05 


300 Myr 


1.28 X 10** 


235 


3.05 


B 





0.1 


120 Myr 


2.56 X 10^ 


235 


3.05 


C 





0.3 


40 Myr 


7.68 X 10^ 


235 


3.05 


Ar 


1 


0.05 


350 Myr 


1.28 X 10^ 


235 


3.05 


Br 


1 


0.1 


200 Myr 


2.56 X 10*5 


235 


3.05 


Cr 


1 


0.3 


40 Myr 


7.68 X 106 


235 


3.05 


B/2 





0.1 


120 Myr 


1.28 X 10« 


58.8 


6.09 


Br/2 


1 


0.1 


120 Myr 


1.28 X 10^ 


58.8 


6.09 


B/8 





0.1 


120 Myr 


3.20 X lO'' 


3.67 


24.4 


Br/8 


1 


0.1 


120 Myr 


3.20 X 10^ 


3.67 


24.4 



8 

UV radiation 

4 t 




ft 1 1 1 1 

UV radiation 



Fig. 1. — Schematic view of the assumed configuration. Rotating gas disk is embedded in a dark halo potential, which is irradiated by the 
ultraviolet background radiation. The light rays are assumed to be perpendicular to the disk. 
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Fig. 2. — Q{= Csk/ttGE) value is plotted as functions of radial coordinate in the disk. In case Q > 1 is satisfied in a certain region of a 
disk, the region is gravitationally stable, otherwise not. Three curves correspond to the three models A,B and C assuming T = T^i^. 
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Fig. 4. — Phase diagram for models A,B and C. Each dot represents each SPH particle. The vertical dotted lines represent njj 
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Fig. 5. — Same as Fig|3] except that the models are Ar, Br and Cr. 
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Model Ar(300Mvr) 



Model Br(I20Mvr) 



Modeler [4aMvr) 




-0.1-0 J il2 0-4 -0.4^.2 i] 0.: 0.4 -0.4 -OJ ^&.2 0.4 
.i'kpc" i'kpc" 2pcpd 

Fig. 7. — Various quantities in the disks are plotted as functions of vertical coordinate, z. The horizontal axes denote 2, whereas vertical 
axes denote density (top), temperature (middle) and HI fraction ym (bottom). These are the snapshots taken at t = 300Myr for model A, 
t = 120Myr for model B and t = 40Myr for model C, respectively. 



1 

0.1 

^ O.I 

10-' 





Cr 




10' 



10' 



10=' 



10^ 



T[K] 



10' 



Fig. 8. — Temperature Probability Distribution Functions(PDFs) are shown for six models Ar,Br,Cr,A,B and C. Horizontal axes denote 
the temperature, whereas the vertical axes show the differential mass fraction AM/Mjotal found within a logarithmic temperature bin 
A(logr) = 0.06. Here M-^^^^i is the total gas mass of the disk. The histograms denote temperature PDFs at t = SOOMyr for models Ar and 
A, t = 120Myr for models Br and B, t = 40Myr for models Cr and C, respectively. 
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Fig. 9. — Density Probability distribution functions(PDFs) are shown for models A,B,C,Ar,Br and Cr. Horizontal axes denote the density, 
whereas the vertical axes show the differential mass fraction AM/Mtot^l found within a logarithmic density bin A(lognH) = 0.18. The 
histogram drawn by thin solid lines denote the initial PDFs in the models. The thick solid lines denote the evolved PDFs at t = SOOMyr for 
models A and Ar, t = 120Myr for models B and Br, t = 40Myr for models C and Cr, respectively. The vertical dotted thin lines show the 
resolution limits of the runs. 




Fig. 10. — PDFs in models B and Br are compared to the results from the low resolution runs(B/2, Br/2, B/8, Br/8). The axes are the 
same as Fig[9] Top three panels :the PDFs for models Br(black histogram) , Br/2(red), Br/8{blue) at t = 40Myr(left), lOOMyr(middle) and 
120Myr(right). Density resolution limits, nn i-csi corresponding to the runs Br, Br/2 and Br/8 are shown as the vertical dotted lines with colors 
same as the histgrams. Bottom panels : PDFs for models B, B/2 and B/8 at three epochs. The colors are same as the top panels. 



